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Abstract 

The bounds on the neutrino mixing angles and CP Dirac phase for an 50(10) model with lop- 
sided mass matrices, arising from the presence of 16// and 16// Higgs representations, are studied 
by variation of the one real and three unknown complex input parameters for the right-handed 
Major ana neutrino mass matrix. The scatter plots obtained favor nearly maximal atmospheric 
neutrino mixing, while the reactor neutrino mixing lies in the range 10~^ sin^ ^13 < 1 x 10^^ 
with values greater than 10~^ most densely populated. A rather compelling scenario within the 
model follows, if we restrict the three unknown complex parameters to their imaginary axes and 
set two of them equal. We then find the scatter plots are reduced to narrow bands, as the mixing 
angles and CP phase become highly correlated and predictive. The bounds on the mixing angles 
and phase then become 0.45 < sin^ 6^23 ^ 0.55, 0.38 < tan^ 6*12 ^ 0.50, 0.002 < sin^ 6*13 < 0.003, 
and 60° ±(5cp ^ 85°. Moreover, successful leptogenesis and subsequent baryogenesis are also 
obtained, with r]B increasing from (2.7 to 6.3) x 10^^^ as sin^ ^23 increases from 0.45 to 0.55. 
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I. INTRODUCTION 

Since the discovery of atmospheric neutrino oscillations Q], more recently that of solar 
neutrino oscillations 21, and still more recently confirmation of these two types of neutrino 

nn 

oscillations with accelerator- and reactor-produced neutrinos p|, |J| , many models have been 
proposed in the literature to explain the mass and mixing parameters associated with these 
oscillations jsj. (We have cited just the most recent references for those experiments and 
several review articles for the models proposed.) Some models restrict their scope to the 
lepton sector, while others such as grand unified models are more ambitious and attempt to 
explain the masses and mixings in both the lepton and quark sectors. The textures for the 
mass matrices obtained in the models may be simply postulated at the outset, derived from 
the observed mixing matrix with a diagonal charged lepton mass matrix, obtained from a 
certain unification group in conjunction with a well-defined family symmetry, or rely solely 
on the Clebsch-Gordon coefficients of that group. 

Until now, all viable models for the lepton sector have only had to explain the apparent 
near maximal ^23 mixing angle for the atmospheric neutrinos and the less than maximal 612 
large mixing angle (LMA) solution for the solar neutrinos, while satisfying the two observed 
mass squared differences, Am^g and Am2i, and the upper bound for the still unobserved 
6'i3 mixing angle. On the other hand, the neutrino mass hierarchy (normal or inverted), 
the Majorana vs Dirac nature of the neutrinos, the appropriate values for ^13 and the CP- 
violating phases (one Dirac and two Majorana) remain undetermined experimentally. 

With the next more precise round of neutrino oscillation experiments which will restrict 
some of the unknowns, many of the more quantitative models proposed are not expected 
to survive the more stringent tests of their predictions. In particular, determination of the 
neutrino mass hierarchy will eliminate models based on an approximately conserved L^, — 
— Lt- lepton number if the hierarchy is observed to be normal, while the conventional 
vpe I seesaw models will be strongly disfavored if the hierarchy is observed to be inverted 
2|. Moreover, small as opposed to very small predictions for ^13 pose a critical test for many 
models. It is thus of utmost importance to sharpen the predictions for each model, so that 
the list of presently viable models can be winnowed down as much as possible when the new 
data is obtained. 

One particularly interesting class of models is that based on the 5*0(10) grand unification 
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group. There are two kinds of minimal models in this class: those based on Higgs repre- 



sentations with dimensions 10, 126, 126, and^ossibly also 120 and/or 210 [8]; and those 
based on 10, 16, 16, and 45 representations 9]. The former type generally has symmet- 
ric and/or antisymmetric matrix elements, while the latter type generally involve lopsided 
matrices for the down quarks and charged leptons. The matrix textures are derived either 
from the Clebsch-Gordon coefficients required or from some Abelian or non-Abelian family 
symmetry. These two types of 50(10) models tend to predict 6*13 angles either near the 
present upper CHOOZ bound of 12°, or in the 1° — 3° range, respectively. 

The purpose of this paper is to sharpen the allowed values for the mixing parameters of 
the lopsided 5*0(10) model proposed some time ago by Babu, Barr, and the author jlol |. 
Upon further refinement by the latter two authors jul, the model has become very predictive 
with all four quark mixing parameters and nine quark and charged lepton masses determined 
by just eight of the twelve model input parameters. The four remaining input parameters 
characterize the right-handed Majorana mass matrix with three of them complex in general. 
In most of the previous papers based on this model, the authors took the latter parameters 
to be real and adjusted them to give good fits to the then known neutrino mass and mixing 
data. Here we allow all three to be complex and obtain scatter plots of the predicted neutrino 
mixing parameters. In so doing, we obtain some bounds on ^13 and 62^ which, if violated, 
will rule out the model which has been highly successful todate. By introducing two very 
small additional parameters in the Dirac neutrino mass matrix, we find that the observed 
baryon asymmetry can be obtained through resonant leptogenesis in the model. In order to 
obtain the latter predictions, the bounds obtained on the neutrino mixing parameters and 
CP phase are even more tightly constrained. With maximal CP violation in the Majorana 
mass matrix and two of the three parameters equal, the mixing angles and phase become 
highly correlated with each other and the baryon asymmetry generated in the model. 

In Sec. II we briefiy review some of the model details and present the results of the 
parameter searches in Sec. III. A summary of the constraints found in the model is presented 
in Sec. IV, where we also briefiy comment on the predictions of the model found by Jankowski 
and Marbury for the /i — > e + 7 lepton fiavor- changing branching ratio. Their prediction 
is also crucial for the viability of the model and has a chance of being tested by the MEG 
collaboration even before the values of the mixing angles are further restricted by new 
reactor and/or long baseline experiments. 
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II. S'O(IO) MODEL AND ITS ARBITRARY PARAMETERS 



The 5*0(10) grand unified model as initially proposed by Babu, Barr and the author 
[lol |. and later updated by Barr and the author to apply to the large mixing angle (LMA) 
neutrino oscillation solution has already been described in great detail in the literature 
Here we just summarize the salient features and refer the reader to those references for more 
information. 

The model is based on the minimum set of Hig gs f ields which solves the doublet-triplet 
splitting problem as proposed by Barr and Raby This requires just one 45h whose 

vacuum expectation value (VEV) points in the B — L direction, two pairs of 16^^, 16jj's 
which stabilize the solution, along with several Higgs fields in the IOh representations and 
Higgs singlets. The Higgs superpotential exhibits the U{1) x Z2 x Z2 symmetry which 
is used for the flavor symmetry of the grand unified model. The combination of VEVs, 
{45h)b~l, (l(16/f)) and (1(T6h)) break 50(10) to the Standard Model, where the latter 
two VEVs point in the SU{5) singlet direction. The electroweak VEVs arise from doublets 
in the 5 and 5 representations of SU{5) in the combinations 



(1) 



{H^) = (5(10^,)), 
{Hd) = (5(10j^))cos7+(5(16'^))sin7, 
while the mixing orthogonal to Ha gets massive at the grand unified theory (GUT) scale. 



Matter superfields appear in three chiral 16's, two pairs of vectorlike 16 and 16, two lO's 
and several I's. The mass matrices follow from Froggatt-Nielsen diagrams in which 
the superheavy fields are integrated out, with the flavor symmetry allowing only certain 
particular contributions to their mass matrix elements. The Dirac mass matrices for the up 
quarks, down quarks, neutrinos and charged leptons are found to be 
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in the convention where the left-handed fields label the rows and the left-handed conjugate 



fields label the columns. 

The entries in the first row and first column are the most uncertain, especially for the up 
quark and Dirac neutrino mass matrices, for they arise from higher order terms involving 
several integrations out of the massive neutrino lines. As originally proposed, the parameters 
and S'j^ were absent. In a later application of the model to leptogenesis [l6|, we have 
observed that the baryon asymmetry lies closer to the preferred value of tjb — 6.2 x 10~^° 



In 



|l7| . if we allow 6n and 5^ to take on nonzero values in the Dirac neutrino matrix M^, 
which are small enough in Mjj, however, so as not to upset the good results in the quark 
sector. We shall consider both possibilities for these two parameters in the searches to be 
presented. 

The lopsided nature of the down quark and charged lepton mass matrices arises from the 
dominance of the a 16^ contribution to the electroweak symmetry breaking over that of 
the 10^. The pronounced lopsidedness with tan/9 ~ 5 readily explains the small Vcb quark 
mixing and near maximal atmospheric neutrino mixing for any reasonable right-handed 
Majorana neutrino mass matrix, M^. All nine quark and charged lepton masses, plus the 
three CKM angles and CP phase, are well-fitted with the original eight input parameters 
defined at the GUT scale to fit the low scale observables after evolution downward from 
Agut ^ 2 X 10^*5 GeV ^|: 

mu — 113 GeV, ^ 1 GeV, 

a = 1.83, e = 0.147, 

(3) 

S = 0.00946, 6' = 0.00827, 

4> = 119.4°, T] =6x 10-^ 

The effective li ght neutrino mass matrix, M^, is obtained from the conventional type I 
seesaw mechanism |l9| once the right-handed Majorana mass matrix, Mr, is specified. While 
the large atmospheric neutrino mixing z/^ ^ Ur arises primarily from the structure of the 
charged lepton mass matrix Ml, the structure of the right-handed Majorana mass matrix Mr 
determines the type of ^ u^^, solar neutrino mixing, so that the solar and atmospheric 
mixings appear to be essentially decoupled in the model. The LMA solar neutrino solution 
is obtained with a special form of Mr yjj], which can be explained by the structure of the 



5 



Froggatt-Nielsen diagrams: 



/ crrf —her] arj \ 
Mr= -ber] e^ -e Ar, (4) 
\ arj —e 1 J 

where the parameters e and r] are those introduced in Eq. Q for the Dirac sector. The 
vanishing 2 — 3 sub determinant of Mr arises naturally if the Yukawa couplings are universal 
at the GUT scale. The new parameters a, b and c are undetermined but expected to be of 
(9(1), since the first row and first column of Mr have been properly scaled by powers of the 
very small rj parameter. As such, the right-handed Majorana mass matrix exhibits a very 
strong hierarchy. 



Given the right-handed Majorana mass matrix above and S 
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Note that the very strong hierarchy in Mr can nearly balance the large hierarchy in M^ 
to yield a rather mild mass hierarchy for the My light neutrino mass matrix. In fact, 
if Mr OC MJjMn were to hold, which has a texture similar to that in Eq. (jU, the seesaw 
formula would yield a light neutrino mass matrix exactly proportional to the identity matrix 
resulting in a three-fold mass degeneracy. The author has argued in Ref. j?! that a normal 
hierarchy solution is much more stable than an inverted hierarchy solution for such type I 
seesaw models, and our results apply for this kind of solution. Barr has recently proposed 
more general possibilities which result in such a mild neutrino mass hierarchy 20|. 

The parameter Kr can be chosen real and to a large degree sets the scale for the atmo- 
spheric neutrino mixing mass squared difference. In most previous studies of the model, 
aside from that of leptogenesis, we chose the remaining three parameters, a, b and c, to be 
real. In order to determine the range of predictions for the neutrino mixing angles, we now 
allow all three of them to be complex. 
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III. RESULTS OF THE PARAMETER SEARCH 



In order to carry out the search of the parameter space for a stable normal hierarchy 
solution, we shall require that the oscillation parameters determined by the model lie in the 
following presently allowed ranges at the 90% confidence level -Q, Q: 

^ml^ = (1.9 - 3.0) X 10-3 eV^, sin^ 623 = 0.36 - 0.64, 

Amj^ = (7.4 - 8.5) x 10"^ eV^, tan^ 6^ = 0.33 - 0.50, 

(6) 

Am^i~Am^2, sin^ 2^13 < 0.16. 

The allowed variation in sin^ 6'23 corresponds to the present bound of sin^ 2^^23 > 0.92 Q|. 
We find that the following bounds for the magnitudes of the input parameters in cover 
the above experimentally allowed ranges: 

|a| < 2.4, \b\ < 3.6, |c| < 3.6, (7) 

where we have set Af> = 2.85 x 10^"^ GeV as explained below. We then use Monte Carlo 
techniques to throw points in the complex a, b and c planes which uniformly cover the disks 
whose radii are indicated above. 

It should be noted that in this paper, we shall neglect any running of the neutrino mass 
and mixing parameters from the GUT scale down to the low scales. This has become more 
or less standard procedure in the literature and is less susceptible to corrections when tan f3 
is low, the neutrino mass hierarchy is normal, and the two low-lying neutrino mass eigen- 



states have nearly opposite CP-parities j2^. All three conditions hold in the model under 
consideration. In fact, we can use the convenient add-on program package for Mathematica 
called "REAP" provided by Antusch et al. j^] to test the assumption that the evolution 
effects are negligible for our case. 

For this purpose, what is required in addition to the mass matrices in Eq. (j2I) is specifi- 
cation of the appropriate Higgs VEVs, Vu = (5(10//)) and Vd = (5(10//) in Eq. (jT)), together 
with the two input scaling masses, rriu and mj:,; from these the corresponding Yukawa cou- 
pling matrices can be determined. Since (Hy) receives a contribution only from (5(10//)), 
the 33 element of the Dirac matrix Mn is given by 

(M^)33 = mu = {Yn)33Vu, (8) 



where Yat is the Yukawa couphng matrix for the Dirac neutrino mass matrix and Vu = v sin (3 
with V = 174 GeV as usual. On the other hand, (Hd) has two contributions from (5(10//)) 
and (5(16i/)), but only the first contributes to the 33 element of the charged lepton mass 
matrix Mi; hence 

{Ml)33 =fnD = (Yl) 33V d cos 'J, (9) 

where Yl is the Yukawa coupling matrix for the charged leptons and cos 7 is the projection 
of onto the VEV of (5(10//)) = Vd = f cos/5. With tan/5 = Vu/vd = 5 as suggested 
earher to provide a sufficiently lopsided charged lepton mass matrix, and the values for mu 
and m/j in Eq. Q, we find that 

Yn = 0.662(Mjv/m//), 

(10) 

Yl = FiCOS7 = 0.0293(Mi/m/)), 

where we have defined Y^ as the effective charged lepton Yukawa coupling matrix due to 
the cos 7 suppression. We assume that the VEVs, f„ and Vd, are held constant while Y^ 
and Y^ evolve from the GUT scale to the weak scale according to the renormalization group 
equations. 

By inserting this information into the REAP package of Antusch et al. Q| , we can obtain 
the evolved results for the neutrino mixing parameters. We find for various test cases that 
the evolved mixing angles and phases are changed by less than 1% from their GUT scale 
values, while the light neutrino masses are reduced by 20%. The latter discrepancy can be 
corrected if we simply reduce the right-handed Majorana scaling factor A/j by 20% from 
2.85 X 10^^ GeV to 2.25 x 10^^ GeV, with essentially no alteration in the mixing angles and 
phases. With this understanding of the role played by evolution of the Yukawa couplings, 
we simply ignore the evolution and proceed to use the GUT scale parameters given earlier. 

Scatter plots are shown in Fig. 1 for the original choice of 77 = 0.6 x 10~^ and Sn = 6'jsf = 0. 
The atmospheric and solar neutrino mixing regions are shown in Figs. 1(a) and (b), where 
Am32 and Am^i are plotted against sin^ 29atm = 4|?7^3p|f/r3p and sin"^ 29soi = 4|f/eip|f/e2p, 
respectively, the atmospheric and solar mixing angles in terms of the neutrino mixing matrix 
elements. For the atmospheric neutrino mixing, we observe from Fig. 1(a) that while An 
sets the scale for the mass squared difference, there is a spread of points for Amig in the 
allowed interval (1.9, 3.0) x 10^^ gY2. gq. (5). A value oi Ar = 2.85 x 10^^ GeV has been 
chosen, in order to center the allowed region predicted by the model on the experimentally 
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allowed region as determined by the Super Kamiokande data. For their best fit point of 
Am|2 = 2.4 X 10^'^ eV^ and sin^ 2^^23 = 1.00, the model prefers a nearly maximal mixing 
angle of sin^ 2^23 ~ 0.98, in agreement with their finding. Note that the highest concentration 
of model points also occurs in this region surrounding their best fit point. For the solar 
neutrino mixing depicted in Fig. 1(b), we observe a nearly uniform spread in Amg^ over 
the region allowed in Eq. (jH)) by the present data. On the other hand, for the solar mixing 
angle there is a slightly greater concentration of points toward the smaller allowed values 
for sin^ '2.0 soi- 

Although it has been customary to present the neutrino mixing data in /S.m? vs. sin^ 29 
mixing planes as in Figs. 1(a) and (b), some authors 2J| have suggested that sin^ ^ is a more 
sensitive measure of the mixing angle, especially for angles close to a maximal mixing of 45°. 
Of particular interest is whether ^23 lies below 45° or above 45°, if the mixing is not exactly 
maximal. With a normal mass hierarchy as suggested by our model, the former signifies that 
the highest neutrino mass eigenstate corresponds more to a v^- flavor composition rather than 
to a composition, while the opposite is true for the latter possibility. In Figs. l(c)-(f) we 
also present results for (c) tan^6'i2, (d) sin^6'i3, (e) (5cp, and (f) the baryon asymmetry, 77^, 
vs. sin^ ^^23- We see from these plots that there is just a slight preference for values of sin^ 62^ 
greater than 0.5 than for those less than 0.5. Concerning the solar neutrino mixing angle 
distribution in Fig. 1(c), the full range of 0.33 < tan^ 6'i2 < 0.50 is populated, but values 
below 0.4 are slightly preferred to those above 0.4 in keeping with remarks made about Fig. 
1(b). Nearly all the allowed points for sin^ ^13 in Fig. 1(d) are found to lie below 0.012 and 
above 1 x 10^^, with those values above 0.001 most densely populated. Interestingly enough, 
the present best fits for sin^ ^13 are 0.006 and 0.004 from three-neutrino mixing analyses with 
all oscillation data considered 2^ • The model seems to be rather insensitive to the leptonic 
CP- violating phase 5c p in this more general search of the complex values for a, h and c, for 
the full range from —90° to 90° is essentially uniformly populated in Fig. 1(e). Finally we 
note that values for the baryon asymmetry in Fig. 1(f) are concentrated in the 10^^'^ — 10~^^ 
range, which are much too small to explain the observed baryon asymmetry in the universe. 

In Fig. 2, the corresponding scatter plots are presented for the choice t] = 1.1 x 10~^, 5^ = 
— l.Ox 10~^ and 5^ = —1.5x10"^. This latter choice has been made as a result of a protracted 
search, so as to maximize the calculated baryon asymmetry as we shall see, without spoiling 
the results in the quark sector. In the previous study reported in the last paper of reference 
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[16[, the corresponding values taken were t] = 0.6 x 10~^, Sn = —0.65 x 10~^ and = 
— 1.0 X 10^^. But without altering the parameter, rj, we observed a maximum of rjB which 
fell short of the presently observed baryon asymmetry. 

The atmospheric neutrino mixing plane in Fig. 2(a) exhibits a somewhat stronger prefer- 
ence for lower values of Amf2 than in Fig. 1(a), while the solar mixing angle distributions in 
Figs. 2(b) and (c) exhibit a slightly greater preference for smaller solar mixing angles. From 
Fig. 2(d) we see the upper bound on sin^ 6'i3 is predicted to lie lower than before in Fig. 
1(d) and is now given by sin^ ^13 < 0.007. The distribution is also somewhat more diffuse 
and less concentrated near the upper limit. Again no preferred values of the CP- violating 
phase 6cp are observed in Fig. 2(e), but the baryon asymmetry is shifted to higher values 
in the broader range, 3 x 10^^^ to 5 x 10^^^ in Fig. 2(f). 

The scatter plots presented in Figs. 1 and 2 were obtained by throwing 400,000 sets of 
points for the parameters a, b and c in the ranges given in Eq. (|7j) and demanding that 
the calculated mixing parameters lie in the experimental ranges quoted in Eq. (jH)). But the 
question arises whether more sparsely scattered points can be obtained outside the apparent 
boundaries. This is of special interest for the baryon asymmetry, where in fact one isolated 
point in Fig. 2(f) does show up with rjB = 2.6 x 10^^°, well outside the major clustering 
but still below the observed value of 6.2 x 10^^'^. Rather than throw 10 - 100 times as many 
points, we searched for input parameters which favor larger values oi rjB- We found 
such a region defined by constraining the parameters a, h and c to lie near or along their 
imaginary axes in the complex planes, still bounded by the magnitudes given in Eq. 0. 

In Figs. 3(a) and (b), we present the scatter plots in the atmospheric and solar neutrino 
mixing planes obtained by throwing 100,000 sets of points for positive and negative imaginary 
values of a, h and c bounded by the magnitudes of Eq. ((7j), again for the new values of r/, 5n 
and 6'^. While the solar mixing plane of Fig. 3(b) is very similar to that of Fig. 2(b), the 
atmospheric mixing plane in Fig. 3(a) is dramatically affected as the spread of points now 
lie along the arc shown. The upper branch corresponds to values of sin^ ^23 < 0.5 and the 
lower branch to values of sin^ ^23 > 0.5. The small secondary upper branch is obtained when 
|c| < |a|, while the major upper branch is found when |c| > \a\. Further insight into the 
restricted results are obtained from Figs. 3(c)-(f). Bands develop in all four plots of (c) 
tan^6'i2, (d) sin^6'i3, (e) 5cp, and (f) rjB vs. sin^6'23. One of the most striking results is 
that sin^ 6*13 is now limited to the restricted range (0.0018, 0.0035). As seen in Fig. 3(f), 



10 



29 points have appeared in the rjB > 1-0 x 10~^" range. 

In Fig. 4 similar plots for the mixing planes and distributions are shown, where 200,000 
sets of points for imaginary values of a, b and c were thrown with the restriction rjB > 
1.0 X 10"^'^ applied. In all 84 points survive this cut. We now see that the upper branch in the 
atmospheric neutrino mixing plane of Fig. 4(a) is essentially eliminated, for sin^^^23 ~ 0.47 
is required for most points to obtain such large values of rjB- The scatter points in the solar 
neutrino mixing plane is little affected aside from the density of points as seen in Fig. 4(b). 
The bands in Figs. 3(c), (d) and (e) become extremely narrow in Figs. 4(c), (d), and (e) 
when the cut r]B > 1.0 x 10""^° is imposed. A strong correlation arises between the solar 
mixing angle in tan^^^i2 and the atmospheric mixing angle in sin^^23; the larger the former, 
the smaller the latter. The upper bands in Figs. 4(d) and (e) are obtained when both 
parameters a and b have positive imaginary values and yield a positive Dirac phase Sep, 
while the lower bands arise from negative imaginary values for a and b. Note that for the 
restriction rjB > 1.0 x 10~^°, the leptonic CP phase lies in the two ranges, ±(50°, 80°). 

The distribution of points in tjb and sin^^23 which survive the cut of tjb > 1.0 x 10~^° 
is shown in Fig. 4(f). Several points with rjB — 6.0 x 10~^° are observed to be clustered 
around sin^ ^23 — 0.6. For these few points, values of c ~ ±a seem to be preferred. This 
suggested that we construct scatter plots with all three parameters, a, b and c taken to 
be imaginary with the restriction c — a. The corresponding plots are shown in Fig. 5, 
where 100,000 sets of points are thrown and 381 points survive. It is rather remarkable that 
the atmospheric neutrino mixing plane in Fig. 5(a) is restricted to a short arc centered at 
Am^2 = 2.45 X 10"^ eV^ and with sin2 2^atm lying in the range (0.985, 0.995), as 0.45 < 
sin^ ^23 < 0.55 with 2.30 x 10"^ < Amj^ < 2.65 x IQ-'^ cV^ for our choice of Ar. The solar 
mixing plane in Fig. 5(b) is also somewhat truncated at the lower end of the spectrum as 
now sin^2^soi ~ 0.80. The bands in Figs. 5(c)-(e) are similar to those in Fig. 4(c)-(e) but 
again somewhat truncated by the maximum value of ^23- The relevant values of sin^ ^13 now 
lie in the interval (0.0020, 0.0031). The CP phase is bounded by ±(60°, 86°) in Fig. 5(f). 

With the strong correlation of the mixing angles noted above, rjB is now directly related 
to the mixing angle ^^23 in Fig- 5(f), as the random scattering in Fig. 4(f) has disappeared. 
The observed baryon asymmetry value of 6.2 x 10"^'^ is obtained for sin^ ^23 — 0.55, which 
implies tan^^i2 c:=^ 0.38, sin^ ^13 ~ 0.0021 and 5cp — ±60°. For this case the sum of Ou 
and the Cabibbo angle is given by 9c + ^12 = 44.8°, very close to the value of 45° which 
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is referred to as quark-lepton complementarity in the literature On the other hand, 

maximal atmospheric neutrino mixing with sin^6'23 = 0.5 corresponds to tan^ ^^12 = 0.44, 
sin^^is ~ 0.0024, Scp = ±75°, and r/^ = 3.8 x 10-^°. 

As an example of a complete solution, we present the results for a special case which arises 
for the choice of right-handed Majorana neutrino parameters: a = 0.5828i, b = 1.76702, c = 
0.5828i with Ar = 2.85 x 10^^ GeV, = 1.1 x 10-^ 6n = -1.0 x 10-^ and 5^ = -1.5 x 10"^ 
One obtains rjB = 6.2 x 10^^° with the following mixing parameters 

Aml^ = 2.30 X 10-3 eV^ sin^ ^23 = 0.550, 

Amj-^ = 7.71 X 10-5 eV^, tan^ ^12 = 0.388, 

sin^^is = 0.0022, 6cp =63°, (H) 

Xi = -55.3°, X2 = 31.8°, 



where Xi ^ind X2 are the Majorana phases in the Majorana phase matrix, $ = 
diag(expixi, expix2, 1)- For this case, ms = 48.9 meV, m2 = 9.30 meV, mi = 3.05 
meV, M3 = 2.91 x 10^^ GeV, Mi ~ M2 ~ 5.40 x 10^ GeV, M2 - Mi = 1.06 x 10^ GeV, and 
Fi = 2.07 X 10'^ GeV where Fi is the width of the lightest right-handed state. 

Further study reveals that the narrow band structures in Fig. 5, obtained with a, b 
and c pure imaginary and c = a, arise because b is limited to the very narrow ranges, 
b = ±(1.7 to 2.1)i, for all 381 points which satisfy the mixing parameter constraints imposed 
in Eq. On the other hand, the ratio a/b varies in the range a/b = (0.33 to 0.53) along 
each band as sin^^^23 decreases from 0.55 to 0.45. The pair of slightly displaced bands in 
Fig. 5(d) is associated with the choices of positive or negative imaginary values for a and b. 
Clearly the imposition of these restrictions for the a, b, c parameter set has greatly enhanced 
the probability that the baryon asymmetry will lie in or near to its observed value. This 
makes for a rather compelling scenario within the model. 

The question arises as to how likely is it that such restricted values should apply for 
the Majorana mass parameters. It is easy to see from the Froggatt-Nielsen diagrams for the 
Majorana mass matrix that if only one Higgs VEV is responsible for breaking lepton number, 
by the U{1) x Z2 x Z2 family symmetry of the model only one diagram will contribute to 
each matrix element, the structure of Mp will be purely geometrical and its determinant 
will vanish, cf. Ref. jl^. With b ^ c = a, only one diagram contributes to each of the 11, 
13, and 31 matrix elements, while two or more diagrams will be required for the 12 and 21 
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elements. This is the simplest possibility for a nonvanishing determinant, as required for 
the seesaw mechanism. The fact that all three parameters are purely imaginary makes for a 
maximum CP violation in the right-handed Majorana mass matrix as all other parameters 
are real. 

To understand the apparent difficulty in reaching values in the observed range oi rjB — 
6.2 X 10~^°, we note the following. The scenario of resonant leptogenesis j27| arises naturally 
in the model, for the two heavy right-handed Majorana neutrinos have nearly equal masses 
and nearly opposite CP-parity. However, one must impose the condition that the mass 
difference of these two heavy Majorana states be at least as large as the half-width of either 
state for the resonance formula to make sense. This forced us to raise the magnitudes of 
?7, 5n and 5'^ in the Dirac neutrino mass matrix. For smaller values of these parameters 
than those chosen in the relevant figures, fewer points reach the 10~^° range of rjB-, as the 
range of sin^ ^23 is shifted lower with more points occuring below 0.50 than above 0.50. It 
is interesting that in going from the set of matrices with the original choice of 77, 5n and 
in Fig. 1 to the new choice of values in Figs. 2-5, the Vus mixing parameter in the quark 
mixing matrix is raised from 0.2220 to 0.2240, in better keeping with the new observed value 
which is preferred for unitarity reasons. The one negative effect in making this change is 
that the up quark mass is raised from 3.5 MeV to 5.5 MeV after evolution from the GUT 
scale, which is above the preferred range. 

While the condition for thermal resonant leptogenesis is easily satisfied for such right- 
handed Majorana masses as in the example given above, the overproduction of gravitinos in 
the early universe is a problem in the supergravity scenario of SUSY breaking; this can be 



28| . As illustrated in 



alleviated if the SUSY breaking occurs via the gauge-mediated scenario 
Ref. Q, this problem can be completely avoided and satisfactory leptogenesis and baryon 
asymmetry obtained, if the model is expanded to include three additional intermediate 
mass scale singlets which lead to a double seesaw mechanism. However, with this expanded 
seesaw mechanism many more parameters are introduced, and the model becomes much less 
predictive. 
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IV. CONCLUDING REMARKS 



From the above results we can conclude that the model with a normal mass hierarchy can 
easily fit the presently observed neutrino mixing data. As originally proposed, eight of the 
twelve model parameters have been fixed by the quark and lepton masses and quark mixing 
data, while the remaining four (three complex and one real) parameters of the right-handed 
Majorana mass matrix are varied. The real A/j parameter centers the Amgj model spectrum 
on the observed 90% confidence level interval given in Eq. (jH)). For the best fit point of 
— 2-4 X 10~^ eV^, the predicted atmospheric neutrino mixing angle is nearly maximal, 
corresponding to sin^ 2^23 ~ 0.98 in agreement with the observed best fit value of 1.00. No 
such preference for a particular value of tan^ 612 in the allowed range is found. The predicted 
values for sin^ 613 are of special interest. There we found a range of 1 x 10^^ < sin^ ^13 < 0.01 
which lies noticeably lower than that predicted by most other models where values very close 
to the upper bound of 0.04 are more generally expected |29|. This is true not only for the 
nearly conserved Le — L^ — Lr models, but also for the 5*0(10) models with symmetric 126// 
and 126h Higgs representations. The lower values predicted in our model have relevance for 
the future study of neutrino oscillations with reactor beams. If Ue disappearance is observed 
very close to the present CHOOZ bound by the next generation reactor experiments, the 
5*0(10) model under consideration will be ruled out. On the other hand, it is interesting to 
note that if future reactor experiments do not see some Uf, disappearance, the model predicts 
it will certainly be seen with superbeam and/or neutrino factory experiments. While there 
remains much interest in determining the Dirac phase Sep for CP violation, it is rather 
surprising that the model is able to accommodate almost any value. 

With the original Dirac neutrino mass parameters, the resonance leptogenesis inherent 
in the model generated a baryon asymmetry which falls more than two orders of magnitude 
short of the observed value. In an effort to increase the effect, the original very small rj 
parameter was adjusted and two new equally small parameters, 6n and 5^, were added to 
the Dirac neutrino mass matrix. This helped to recover two missing orders of magnitude 
for rjB, provided the three complex parameters in the right-handed Majorana mass matrix 
were chosen very near to or on their imaginary axes. By selecting only those points which 
satisfied the observed neutrino mixing results and for which tib ^ 1.0 x 10"^*^, we observed 
that the mixing angles lie along very narrow bands. 
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These bands were further dramatically tightened by choosing not only a, h and c purely 
imaginary, but restricting c to be equal to a. For this situation a limited range of 0.45 < 
sin^^23 ^ 0.55 was found over which tan^ 6*12 falls from 0.50 to 0.38, sin^^is varies from 
0.0031 to 0.0020, and the Dirac CP phase varies from ±(85° to 60°). On the other hand, 
a narrow band also developed for the baryon asymmetry due to the strong mixing angle 
and CP phase correlations, as it rises from (2.7 to 6.3) x 10~^° with increasing sin^ ^23- The 
appearance of such narrow bands is directly related to the very narrow acceptable range 
for 6, while the ratio a/h varies along the length of each band. We argued that the c = a 
restriction can be understood in terms of one Froggatt-Nielsen diagram for each of the 11, 
13, and 31 matrix elements of M/j, while two such diagrams can contribute to the 12 and 
21 elements. All these results are of considerable interest and will severely test the model 
in the future. 

Since the light neutrino mass spectrum has a normal hierarchy, the model predicts no 
neutrinoless double beta decay should be observed in the near future. In fact, for all choices 
of the model parameters which lead to acceptable ranges for the mixing parameters of Eq. 
dni), we find the effective mass is mpp ~ 0.4 x 10~^ eV jp], three orders of magnitude below 
the present bound {3^ . 

In conclusion, we address briefly the issue of lepton flavor violation in the decay process 



— > e + 7. Jankowski and Maybury 



13 have studied this particular class of 5*0(10) models 



with lopsided charged lepton mass matrices in order to determine the expected branching 
ratio for the above decay, under the assumption that the supersymmetric grand unified 
model breaks directly to the constrained minimal supersymmetric standard model. Given 
the combined constraints on the CMSSM parameters from direct searches and from the 



WMAP satellite observations 
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3l|, they find that for a low value of tan/5 ~ 5 preferred 



by the model, the branching ratio should be of (9(10^^^) and very close to the present 
upper limit of 1.2 x 10~^^ [2^. Hence the /i ^ e + 7 decay should be observed in the new 
MEG experiment which is expected to start taking data in 2006. This prediction is 
based on the slepton masses being in the 125 - 150 GeV range, while the neutralino masses 
are around 300 - 500 GeV. These ranges are also favored in CMSSM models b y th e recent 
study of the best chi-squared fits to all the data with tan/5 ~ 10 by Ellis et al. [33[. On the 
other hand, minimal 5*0(10) models leading to symmetric or antisymmetric mass matrix 
elements favor a smaller lepton- violating branching ratio of (9(10~^^) or less [2^. The 
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coming MEG experiment which is expected to be sensitive down to the 10~^ branching 
ratio level may thus be able to rule out one of the two classes of minimal SO{10) models 
before the new reactor experiments can be launched. 

The author thanks Uli Haisch and Enrico Lunghi for their help with the scatter plots. The 
author also thanks the Theory Group at Fermilab for its kind hospitality. Fermilab is op- 
erated by Universities Research Association Inc. under contract No. DE-AC02-76CH03000 
with the Department of Energy. 
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1. Plots of (a) Am32 vs. sin^2^23, (b) Am^^ vs sin^20i2; and (c) t&n^ 6^, (d) sin^^i2, 
(e) Sep, and (f) r]B vs. sin^ 023 for the original choice of = 0.6 x 10^^ and 
Sn = S'r^ = in the Dirac neutrino mass matrix. A value of = 2.85 x 10^^ GeV 
applies for all figures. 
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Fig. 3. Plots of (a) Aml2 vs. sin^ 26*23, (b) Am^i vs sin^20i2; and (c) t&n^ 6^, (d) sin^^i2, 
(e) Sep, and (f) r]B vs. sin^ 023 for the choice of = 1.1 x IQ-^, 6^ = —1-0 x IQ-^, 
and S'f^ = —1.5 x 10~^ in the Dirac neutrino mass matrix. The parameters a, b, and 
c in the right-handed Majorana mass matrix are constrained to he along their 
imaginary axes. 
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